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Abstract 

We consider the asymmetric exclusion process (ASEP) in one dimension 
on sites i = 1, N, in contact at sites i = 1 and i = N with infinite par- 
ticle reservoirs at densities p a and pf,. As p a and pt are varied, the typical 
macroscopic steady state density profile p(x), x G [a, b], obtained in the limit 
N = L(b — a) —* do, exhibits shocks and phase transitions. Here we derive 
an exact asymptotic expression for the probability of observing an arbi- 
trary macroscopic profile p(x): P N ({p(x)}) ~ exp[-LJ r [ 0)6] ({p(a;)});p a ,p 6 ], 
so that J- is the large deviation functional, a quantity similar to the free en- 
ergy of equilibrium systems. We find, as in the symmetric, purely diffusive 
case q = 1 (treated in an earlier work) , that J~ is in general a non-local func- 
tional of p(x). Unlike the symmetric case, however, the asymmetric case 
exhibits ranges of the parameters for which J r ({p(x)}) is not convex and 
others for which J 7 ({p(x)}) has discontinuities in its second derivatives at 
p(x) = p(x); the fluctuations near p{x) are then non-Gaussian and cannot 
be calculated from the large deviation function. 



1 Laboratoire de Physique Statistique, Ecole Normale Superieure, 24 rue Lhomond, 

75005 Paris, France; email derrida@lps.ens.fr. 

2 Department of Mathematics, Rutgers University, New Brunswick, NJ 08903; email 

lebowitz@math.rutgers.edu, speer@math.rutgers.edu. 

3 Also School of Mathematics, Institute for Advanced Study, Princeton, NJ 08540. 
4 Also Department of Physics, Rutgers. 



1 



Key words: Large deviations, asymmetric simple exclusion process, open 
system, stationary nonequilibrium state. 



2 



1 Introduction 



Stationary nonequilibrium states (SNS) maintained by contact with infinite 
thermal reservoirs at the system boundaries are objects of great theoretical 
and practical interest [[[], ||, |||, |], [|, [f| fj], ||. One is tempted to think 
that, at least when the gradients and fluxes induced by the reservoirs are 
small, the full system behavior is just that of a union of subsystems, each 
in local equilibrium, with spatially varying particle and energy densities or 
equivalently local chemical potential and temperature. This is, however, 
not the entire story, as is clear when one considers the paradigm of such 
systems, a fluid in contact with a thermal reservoir at temperature T a at 
the top and one at temperature T& at the bottom, the Rayleigh-Benard 
system ||. In this system there are long range correlations not present 
in equilibrium systems, which have been measured by neutron scattering 
experiments. This system exhibits, when T& — T a exceeds some positive 
critical value, dynamic phase transitions corresponding to the formation of 
different patterns of heat and mass flow as the parameters are varied. These 
are due to macroscopic instabilities, caused by gravity, but are not derivable 
at present, despite various attempts |J, in terms of a microscopic theory 
of SNS such as that provided by statistical mechanics in the case T a = T5, 
when the system is an equilibrium one. 

In this paper we study the SNS of a model system which, despite its 
simplicity, has some phase transitions and exhibits many phenomena, such 
as long range correlations and non-Gaussian fluctuations, which are very 
different from those of systems in local thermal equilibrium. For this system 
the weights of the microscopic configurations are known and the typical 
behavior has been deduced from them (T^j. Here we determine — again from 
the microscopic weights — the probabilities of various atypical macroscopic 
behaviors, that is, of large deviations. 
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The model we consider is the SNS of the open asymmetric simple exclu- 
sion process (ASEP) JlO| , 11]: a lattice gas on a chain of N sites which we 
index by i, 1 < i < N. At any given time t, each site is either occupied by a 
single particle or is empty, and the system evolves according to the following 
dynamics. In the interior of the system (2 < % < N — 1), a particle attempts 
to jump to its right neighboring site with rate 1 and to its left neighboring 
site with rate q (with < q < 1). The jump is completed if the target 
site is empty, otherwise nothing happens. The boundary sites i = 1 and 
i = N are connected to particle reservoirs and their dynamics is modified as 
follows: if site 1 is empty, it becomes occupied at rate a by a particle from 
the left reservoir; if it is occupied, the particle attempts to jump to site 2 
(succeeding if this site is empty) with rate 1. Similarly, if site N is occupied, 
the particle may either jump out of the system (into the right reservoir) at 
rate (5 or to site N — 1 at rate q. 

More generally, one can also consider the ASEP with the above dynam- 
ical rules supplemented by an output rate 7 at % = 1 and an input rate 5 
at i = N. However, the calculations in the case of nonzero 7 or 5 are more 
complicated, and for that reason we limit our analysis in the present paper 
to the case 7 = 5 = 0. 

It is convenient to introduce the two parameters 

a 1 P ft t\ 

Pa = z , Pb = 1 - -, ; (1-1) 

1 — q 1 — q 

we require that < a, (3 < 1 — q, so that < p a , Pb < 1. These parameters 

have a natural interpretation as reservoir densities. In particular, when 

Pa = Pb the steady state of the system is a Bernoulli measure at constant 

density p a , that is, each site is occupied independently with probability 

p a ; this may be seen for example from the so-called "matrix method" (see 

Section |j). In general, then, we interpret p a and pb as the densities of 

particles in the left and right reservoirs, respectively. 
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The goal of the present work is to calculate the large N behavior of 
Pn({p(x}), the probability of seeing a macroscopic profile p(x) for a < x < b 
in the steady state for a system of N = L(b — a) sites. This probability 
Pn({p(x}) can be thought of as the sum of the probabilities of all microscopic 
configurations such that in each box of L dx sites (with dx <C 1 and Ldx^> 
1), the number of particles is close to Lp(x) dx. The ratio log(P/v({/o(x}))/L 
has a well defined limit for large L, 

logP/v({p(x}) _ . 
hm = -^[afldpix)}; Pa, Pb) , (1-2) 

which depends on b — a, on the density profile p(x), and on the reservoir 
densities. T is called the large deviation functional (LDF) of the system. 

When p a = pb, the steady state is a Bernoulli measure, as described 
above. This measure is just the equilibrium state of a system of parti- 
cles, noninteracting except for the hard-core exclusion, at chemical potential 
log(p a /(l — p a ))- For this system the LDF can be computed by elementary 
means, and is given by |l3|, |l4|, [l5[| 



F[a,b]({p(x)Y, Pa, Pa) 



, \ , P( X ) /-, , \M 1 _ P( X ) 

p(x) log ^-L + 1 - p(x)) log - 

Pa I - Pa 



dx 
(1.3) 



The present paper is devoted primarily to the derivation of the exact ex- 
pression of T in the case p a ^ pb- 

1.1 Additivity and large deviations in the ASEP 



In our earlier work 16, 15] on the large deviation functional for the symmet- 
ric simple exclusion process, corresponding to equal jump rates to the left 
and right, i.e., to q = 1, we were able to use the matrix method JT|, |T| ||, ||] 
to calculate directly the probability of a given macroscopic profile p(x) by 
summing the probabilities of all configurations corresponding to that pro- 
file. For the asymmetric model, however, such a direct calculation is a priori 
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more complicated. For that reason, we follow here a different path, which 
has its origin in an a posteriori observation made in [15]. We noted there 



that while the large deviation functional for the symmetric case is nonlocal, 
it possesses a certain "additivity" property. For the ASEP we first derive an 
additivity property, similar to that of the symmetric model, and from that 
obtain T ' ■ The derivations are given in Section ||. 

The addition formula involves a function Tt related to J- by 

W[a,6]({pO»0};Pa,P&) = F[a,b]({p( X )}'> Pa' Pb) + (P ~ a ) K (Pa, Pb) , (1-4) 

where K(p a ,p ) does not depend on p(x). Since the dynamics in the bulk 
is driven from left to right (for < q < 1), the roles played by the left and 
the right reservoirs are not symmetric and the additivity relation and 
the expressions for K(p a , p ) in (|l.4j ) and for J r t a u({p(x)}; p a , Pb) depend on 
whether p a > p b or p a < p b . 

1.1.1 The case p a > p b 

When p a > p the constant in (|1.4|) is 



K(p a ,Pb)= sup log[p(l-p)], (1.5) 

Pb<P<Pa 

and the additivity relation, obtained in Section |5| below, is that for any c 
with a < c < b, 

K[a,b]({p(?)}] Pa, Pb) (1-6) 

= sup [H[a }C ]({p(x)}; pa, Pc) +'H[c,b]({p{x)};pc,Pb)] ■ 

Pb<Pc<Pa 

Equation (|l.6| ) expresses a relation between the large deviation function of 
the whole system and those of two subsystems connected at the break point 
to a reservoir at an appropriate density p c . 

Once we have the additivity relation (|1.6| ), the derivation of the large 
deviation functional is simple, and we give it here. We divide our system 



6 



into n parts of equal length and apply Ql,6| ) n times; this will introduce 
intermediate reservoir densities p a = po > p\ > ■ ■ ■ > p n = p b . For very 
large n, most of the intervals must have reservoir densities pk-i,Pk at their 
boundaries which are nearly equal, and the LDF for these intervals is ap- 
proximately given by ( |1,3| ) (with p a there replaced by pk-i — Pk)- On the 
other hand, the total length of the intervals for which this is not true will 
approach for large n. Now taking the n — > oo limit and introducing a 
function F(x) as the interpolation of the values po,pi, ■ ■ ■ ,p n , we are lead 
directly to a formula for the large deviation functional: 

F[a,b]({p( x )}'i Pa> Pb) = ~(b- a)K(p a ,p b ) (1.7) 
+ sup/ dxp(x)log[p(x)(l- F(x))] + (l- p(x))log[(l- p(x))F(x)\, 

F(x) J a 

where the supremum is over all monotone nonincreasing functions F(x) 
which for a < x < y < b satisfy 

Pa = F(a) > F(x) > F{y) > F(b) = p b . (1.8) 

This supremum is achieved at a certain function F p (x). Note that without 
the constraints ( |1.8[) one would have F(x) = 1 — p(x). The monotonicity 
requirement, however, makes the determination of F p more subtle (see Sec- 
tion [O]) and the expression of T nonlocal. This will be at the origin of most 
of its interesting properties. 

1.1.2 The case p a < p b 

When p a < p the constant in (|l.4[) is 

K{p a ,Pb) = min[log / o a (l - p a ),logp 6 (l - p b )} , (1.9) 
and the additivity relation, again derived in section [j| is 

H[a,6]({p(^)}; Pa, Pb) = 

mm [H [aiC] ({p(x)};p a ,p c ) +H [Ctb] ({p(x)};p c ,p b )] . (1.10) 

pc—pa,Pb 

7 



By an argument similar to that which led to ( |1,7| ), we obtain the formula 

F[a,b]{{p{x)}\ Pa, Pb)) = ~(b- a)K(p a , Pb ) + (1.11) 

inf \[ dxp(x)log [p{x){\ - p a )\ + (1 - p{x))\og [(1 - p(x))p a ] 
a<y<b {J a 

+ / dxp(x)log [p(x)(l - p b )\ + (1 - p(x))log [(1 - p(x))p b ] 
Jy 

The fact that for each p(x) one has to find in ( |l.llj ) the infimum over y 
makes, in this case too, the large deviation function nonlocal. 

It is interesting to note that the formulas ( |1.7| ) and ( jl.ll ) for the large 
deviation function do not depend on q. 



1.2 Outline of the paper 

In Section |2| we review briefly some known results on the open ASEP, em- 
phasizing the phase diagram. In section || we give a summary of our results 
which follow as consequences of the formulas (1.7) and (1.11) for the large 



deviation function. In Section || we recall the matrix method for carrying 
out exact calculations in the ASEP and give several results obtained by this 
method which are relevant to our considerations here. Then in Section [5] we 



give the derivation of the addition formulas (|l.6| ) and fll.lC ) . An unexpected 



consequence of our results, discussed in Section [| is that the correlations in 
the steady state are not related in any simple manner to the large deviation 
function. We also show in that section that the fluctuations of the number of 
particles in any box of size Lx, with < x < 1, are not Gaussian for certain 
ranges of the parameters p a ,Pb- Section [7| gives some concluding remarks, 
and certain more technical questions are discussed in appendices. 
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Figure 1: The phase diagram of the open ASEP 

2 The steady state of the ASEP with open bound- 
ary conditions 

We describe here the phase diagram of the ASEP with open boundary con- 



ditions, which has been obtained by various methods [20, 12, 18, [19], 21]. As 
indicated above, we consider here only the case 7 = <5 = 0, 0<a, (3 < 1 — q. 

The phase diagram is given in Figure 1, where we have chosen as pa- 
rameters the densities p a and pj, ( |1.1| ) of the two reservoirs. There are three 
phases: a low density phase A with a constant density p(x) = p a in the bulk, 
a high density phase B with a density p(x) = pb, and a maximal current 
phase C with a density p(x) = 1/2. The current in each phase is given by 
J = (1 — (/)p(l — p). The transition lines between these phases are second 
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order phase transitions where p(x) is a continuous function of p a and pb, 
except for the boundary S (p a = 1 — pb < 1/2) between phase A and B, 
where the transition is first order: p(x) jumps from p a to pb- On the line 
S the typical configurations are shocks between phase A with density p a at 
the left of the shock and phase B with density pb at the right of the shock: 

Py{x) = p a &(y - X) + p b @{x -y), (2.1) 

with the Heaviside function. The position y of the shock is uniformly 
distributed along the system [22, 23], and as a result the average profile is 
linear: (p(x)) = p a (l - x) + p b x. 

This phase diagram can be understood easily in heuristic terms. First 
consider an infinite one dimensional lattice on which the initial configuration 
is a Bernoulli distribution at density p a to the left of the origin and pb to the 
right of the origin. If p a < pb, this initial condition produces a shock moving 
at velocity (1— p a — pb)(l— q); thus in the long time limit, the distribution near 
the origin is Bernoulli, with density p a if p a +Pb < 1 and density pb if p a +Pb > 
1. On the other hand if p a > pb, the profile becomes a rarefaction fan: 

p(x,t) = p a if X < X a (t), p(x,t) = Pa + {pb~ Pa){x-X a {t))/(x b {t)-X a {t)) if 

x a (t) < x < Xb(t), and p(x, t) = pb if Xb{t) < x, with x a (t) = (1— q) {1—2 p a )t, 
a = a,b. If 1/2 < pb < p a then the entire fan moves away to the left, if 
Pb < Pa < 1/2 then it moves away to the right, and if pb < 1/2 < p a then 
the origin remains in the fan for all time. These three cases give rise in the 
long time limit to Bernoulli distributions at the origin, with densities p a , pb, 
and 1/2, respectively (see e.g. [^4|] and references therein). 

If now we consider the finite system with left and right reservoirs at 
densities p a and pb, and start with a Bernoulli distribution at density p a at 
the left of some point in the bulk far from the boundaries, and pb at the right 
of this point, then the evolution will be the same as in the infinite system 
until the shock or the fan reaches the boundary, leading, for the asymptotic 
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density in the bulk, to what is given in the phase diagram. This idea is at 
the basis of what has been done recently by Popkov and Schiitz [25, 26] to 
predict boundary induced phase diagrams in more general cases. 

The dashed line p a = pb (a + (3 = 1 — q) in Figure 1 separates what we 
will call the shock region p a < pb and the fan region p a > pb- On this line 
the measure reduces to a Bernoulli measure at density p a (see Section H) 



and the large deviation function is given by (1.3). The line plays no role in 
the phase diagram for the typical profile p but separates phases A and B 
into two subphases, Ax,A% and Bi,B2, which, as we have seen in Section [l], 
can be distinguished by the different expressions (1.7), (1.11) for the large 
deviation function in these regions. 

3 Consequences of the large deviation formula for 
the ASEP 

In this section we describe some consequences of formulas ( |1.7[ ) and ( |1.11| ) 
for the large deviation function. It is convenient to write ( |1.5| ) and ( |1.9| ) in 
the unified form 

K(p a ,Pb) = log p(l - p), (3.1) 

where p = p(x) depends on p a and pb and is obtained from the phase dia- 
gram, Figure 1. Note that p is in fact independent of x except on the line 
S; there p can be any shock profile p y {x) ( j2.1|) , but the value of K(p a ,pb) 
is independent of the shock position y since p a + pb = 1 on S. We also 
introduce the notation 

h(r, f;p)=r log - + (1 - r) log -j— - + log (3.2) 



so that ( |1.7|) and ( 1.11 ) become respectively 



F[a,b]({p( X )h Pa, Pb) 
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pb pb 

= sup / dx h(p(x), F(x)\ p) = \ dx h(p(x), F p (x); p), (3.3) 

F(x) J a J a 

for p a > Pb, and 

F[a,b]({p(x)h Pa, Pb) 

inf { f dx h(p(x),p a ;p) + [ dx h(p(x) , p b ; p)\ , (3.4) 



a<y<b [J a Jy j 

for p a < p , where again p is determined from p a ,Pb through the phase 
diagram. Note that h(r, /; p) is strictly convex in r for fixed /, p, with a 
minimum at r = f. 

3.1 Construction of the function F p 

There is a rather simple way of constructing the optimizing function F p (x) 
in ([T^). Let G p (x) be defined for a < x < b by 

G p (x) = Concave Envelope I J (1 — p(y))<iy| ; (3.5) 

then F p is obtained by cutting off G' p (x) at p a and pb'. 

p a , if G'{x) > p a . 
F P (x) = l G'(x), if Pb <G'(x)< Pa , (3.6) 
. Pb, if G'(x) < p b , 
This construction is verified in Appendix 

Suppose, for example, that p{x) is a constant profile: p(x) = r. Then 
the expression in brackets in ( |3.5[ ) is concave and hence equal to G p (x), so 
that G' p (x) = 1 — r and F p is constant: F p = p a if r < 1 — p a , F p = 1 — r if 
1 — p a < r < 1 — pb, and F p = p if 1 — pb < r. In particular, taking r = p 
we find that 

Fp = p. (3.7) 

For (i) in region A 2 , where 1/2 > p a > p b , r = p a < 1 - p a and so F p = p a = 
p; (ii) in region C, where p a > 1/2 > p , r = 1/2 and so F p = 1— r = 1/2 = p; 
(iii) in region B 2 , where p a > Pb > 1/2, r = pb > 1 — Pb and so F p = p b = p. 
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3.2 The most likely profile 

We will show in this section that the most likely profile is always given 
by p(x) = p; specifically, that F[ a ,b]({p( x )}', Pa, Pb) > for all p(x), that 
F[a,b]({pY, Pa, Pb) = 0, and that F[a,b]({p( x )Y, Pa, Pb) > if p(x) / p. This 
is of course expected and in fact represents an alternate way to obtain the 
phase diagram of Figure 1. 

To obtain the most likely profile in the fan region p a > pb, we note that 
(|3,3|), Q3.7| ), and the convexity of h(r, /; p) in r imply that 

F[a,b]({p( x )Y, Pa, Pb) = / dxh(p(x),F p (x);p) 

J a 

> / dxh(p(x),Fp-;p) 



> 



f dxh(Fp,F p ;p) = 0. (3. 

J a 



Moreover, if p(x) = p then equality holds throughout (|3.8;), so that the 
minimum value of T is zero and p is a minimizer; otherwise the second 
inequality is strict, from which it follows that this minimizer is unique. 
In the shock region p a <C p^, we observe that for fixed y : a ^ y ^ 6, the 



right side of (|3.4| ) is minimized by the unique choice p(x) = p a for a < x < y, 
p{x) = pi for y < x < b. Minimizing over y then implies that, except on the 
first order line S, the optimal profile is again constant with value p, and is 
unique; again the corresponding minimum value of T is zero. On S all values 
of y give the value zero for J 7 , so that the shock profiles p y (x) of ( |2.l| ) form 
a one parameter family of minimizing profiles. Note that a knowledge of the 
large deviation functional is not sufficient to determine the distribution of 
the shock position y mentioned in Section 0. 
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3.3 Convexity 

In the fan region p a > p b , the LDF J T [ a ^({p(x)}; p a , Pb) is a strictly convex 



functional of p(x), since by (3.3) it is the maximum (over the functions F(x)) 
of strictly convex functionals of p(x). This is also true in the symmetric case 
]16| , |l~5| ] and in equilibrium systems not at a phase transition. In the shock 
region p a < p b , on the contrary, T is not convex. This is most easily verified 
on the line S, since it follows from Section |3.2| that on S a superposition 
of minimizing profiles fl2.1|), p{x) = \p y (x) + (1 — X)p z (x), y ^ z, satisfies 



J : ({p(x)}) > for < A < 1. But we will also see in Section 3.5 below that 
for every p a , Pb there is a constant profile p(x) = r* near which T is not 
convex. 



3.4 Suppression and enhancement of large deviations 

The LDF in the fan region p a > p b has similarities besides convexity to the 
LDF in the symmetric case. In particular it is easy to see from ( |1.7| ) that 

F[a,b]({p(z)}\ p a , Pb) > ^ b ]({p(%)Y,p), itpa>Pb, (3.9) 
where we define 

r b r n(r\ 1 - nf^rM 

dx ; (3.10) 



^ b] ({p(x)hp) 



p(x)log^ + {l- p(x))log' 



p 1 - p 

this is the LDF for an equilibrium system at density p (see ( |1.3[ )) But in the 
shock region this inequality is reversed: 



•?>,&] ({pO)}; Pa, Pb) < ^ b] ({p(x)};p), if p a < pb, 



(3.11) 



as can be derived from (|3.1[) and (|3.4j) , since in region B\ (^4i), taking y = a 
(y = b) on the right side of (|3,4j) gives ({p(x)}; p)- On the line S, where 
p = p y (x) (see (|2.1[) ) for some y, ( |3. 11 ) holds for all values of y. Physically, 
Q and (pll ) mean that the probability of a macroscopic deviation from 
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the typical density profile is reduced in the fan region, and increased in the 
shock region, compared with the probability of the same deviation in an 
equilibrium system with the same typical profile. 



3.5 Flat density profiles 

One can readily compute the LDF for constant profiles p(x) = r. In the fan 



region p a > pb one finds, as discussed in Section 3.1, that F p is also constant, 



with F p = p a if r < 1 - p a , F p = 1 - r if 1 - p a < r < 1 - p b , and F p 



Pb 



if 1 — pt, < r. There are correspondingly three different expressions for the 
LDF F{r) = -F[a )6 ](r; p a , p h )\ 



(b-a)h(r, p a ;p), 
(b-a)h(r,l - r;p), 
k (b- a)h(r,p b ;p), 



if r < 1 - p a , 

if 1 - Pa < r < 1 
if 1 — pi, < r. 



Pb, 



(3.12) 



If p a > 1/2 > p , so that p = 1/2, (phase C of Figure 1) and 1 — p a < 
r < 1 — pb, then from ( p. 12 ), 



£(r) = 2JP M (r; 1/2, 1/2) = 2^ 6] (r; 1/2), 



(3.13) 



where again (see ( 3.10 )) 3^ab]( r '> V^) i s the large deviation function for ob- 
serving the uniform density r in a Bernoulli measure with density 1/2. In 
particular, for p a = 1, pb = 0, the probability of observing all sites empty, 
r = 0, or all sites occupied, r = 1, is 2~ N for the Bernoulli measure, where 
N = L(b — a) is the number of sites, so that from ( |3.13 ) the corresponding 
probability for the ASEP is given to leading order by 4r N . 

In the shock region p a < pb, the minimizing y in ( |1.11| ) is y = b if r < r* 
and y = a if r > r* , where 



log 



1— Pb Pa 



(3.14) 
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so that 

(b- a)h(r,p a ;p), if r < r*, 
(b- a)h(r,p b ;p), \ir>r*. 

At r = r* the derivative is discontinuous, with 



(3.15) 



dr 



> 



r— >r*— 



fir 



(3.16) 



r— »r*+0 



Thus ^"(r) is not convex near r = r*, and hence J 7 \ a u({p(x)}; p a , p b )) is not 
convex in a neighborhood of p(x) = r* . 



3.6 Distribution of the total number of particles 

The probability Pl(M) that there are a total of M = rN particles in the 
system can be obtained from the LDF as 

1 



L 



P L (rN) ~ F(r) = F [a>b] ({p(x)}; p a , p b ), 



where p{x) is the most likely profile under the constraint 

r-b 

p(x) dx = r(b — a). 



(3.17) 



(3.18) 



It will be shown in Appendix |B| that p(x) is the constant profile p(x) = r, 
and correspondingly J-(r) = F{r), except in that portion of the shock region 
in which r satisfies 1 — p b < r <1 — p a , where 

1- p b , if x < y r , 
I- p a , iix>y r , 

with 



p(x) 



(3.19) 



6(1 - Pa-r)- q(l - Pb-r) 
y r = . (3.20) 



Pb ~ Pa 

Then from (13.41) it follows that for these values of r, 



F(r) = (b- o)[rlog(l - p a )(l - p b ) + (1 - r) log p aPb - K(p a ,p b )]. (3.21) 
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Remarks: (a) Equation ( |3.19| ) is easy to interpret on the first order line 
S in the phase plane, where p a + p& = 1: it corresponds there to a typical 
shock configuration p yr (x), with y r determined by ( |3.18 ). 



(b) Although, as observed in Section 3.5, T is not convex (at least in the 
shock region) , T is always convex, and is in fact the convex envelope of T. 

(c) We will discuss in Section 3.2 how ( |3.21| ) may be derived by a direct 



calculation from the matrix method. 
3.7 Small fluctuations 

It is natural to ask about the connection between the LDF, which gives the 
probabilities of macroscopic deviations from the typical density profile, and 
the distribution of small fluctuations, i.e., those of order 1/y/N. (In what 
follows we will refer to these simply as "fluctuations." ) In the symmetric case 
discussed in (jDJ , as in an equilibrium system not at a phase transition (in any 
dimension, with N being the number of sites in the system), the distribution 
of fluctuations can be obtained from J 7 as a limit. More precisely if we write 
p(x) = p(x) + Aju(x) and then expand T to second order (the first order 
term being zero) we get a Gaussian distribution for u{x) with covariance 
C(x,x'), where C~ 1 (x,x') = 5 2 J r /5p(x)5p(x') evaluated at p = p. This 



covariance is the suitably scaled microscopic truncated pair correlation [15]. 

For the asymmetric case discussed in this paper, however, the distri- 
bution of small fluctuations need no longer be given by the LDF, as we 



discuss in Section 6A. In fact we show there that 5 2 J 7 /5p(x)5p(x') is dis- 
continuous at p = 1/2 in the interior of region C of the phase diagram, i.e., 
where p a > 1/2 > p^. Furthermore, the fluctuations in this region are no 



longer Gaussian; in Section 6.3 we show this by computing explicitly the 
non-Gaussian distribution of the fluctuations of the number of particles in 
a box of size Ly, with < y < 1. 
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4 The matrix method for the ASEP 



The steady state properties of the ASEP with open boundaries can be cal- 
culated exactly in various ways gg, here we describe the so-called 



matrix method 12, [17], 18, 19], which we will use in the derivation of the 
additivity relations in Section ||. Let us consider two operators denoted by 
D and E, a left vector {W\ and a right vector \V), which satisfy the following 
algebraic rules: 

DE - qED = D + E , (4.1) 
0D\V) = \V), (4.2) 
(W\aE = (W\ . (4.3) 

Any matrix element of the form (W|liY2 • • • Yk\ V) I {W\V) , where Yi denotes 
D or E, can be calculated from these rules (without the need of writing down 
an explicit representation). Thus from ( |4,2j]4,3|) , one has (W|.D| V) / (W\ V) = 
1/(3 and (W|i?|V)/(W|V) = 1/a, and the matrix element of any product of 
n matrices can be expressed through ( |4.lD ~( [0| ) in terms of sums of elements 
of shorter products. 

For the open ASEP as described in Section [l], the probability P({Tj}) of 
the microscopic configuration {t^} in the steady state can be written as jl^] 

p( {T ,» = wnLi^.+^i-T,))^ 

where u = 1 or indicates whether site i is occupied or empty and the 
normalization factor Zn is given by 

Z N = (W\(D + E) N \V). (4.5) 

On the dashed line p a = pb of Figure 1, there exists a one dimensional 
representation of the matrix algebra; it then follows immediately from 
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that the invariant measure is Bernoulli on this line. There are other lines in 
the p a , Pb plane where there exist finite dimensional representations [29. pFj] 



of the algebra ( |4.1| )-( fi~3| ), but the large deviation function has no particular 
or remarkable expression along these lines, and so these cases will be treated 
together with the general case. 

From the algebra ( |4.1]) -( [OD all (equal time) steady state properties can 
(in principle) be calculated. For example the average occupation (rj) of site 
i is given by 

_ (W\{D + Er l D{D + E) N -*\V) ^ 
Zn 

and the two point function is, for i < j, 

(W\(D + Ef- l D{D + Ey~ i ~ 1 D(D + E) N ^\V) 

\ T iTj) = - . (4.7) 

An 

The average steady state current Jn is the same across any bond. It has 
the form 



Jn = - r i+1 ) - q(l - Ti)T i+1 ) 

(W\(D + E) i - 1 {DE - qED)(D + E^-^V) Z N -\ 



.(4.8) 



Zn Zn 
The probability QN lt ..N k (Mi, ..Mk) that there are exactly M\ particles on 
the first Ni sites, M2 particles on the next N2 sites, . . . , particles in the 
rightmost Nk sites is given by 

QN u ...,N k {Mi, . . . ,M k ) - , (4.9) 

An 



with 



1 



2tt 



Y p = — I d6 e~ ieM v (De ie + E) Np . (4.10) 
2iir J 

Clearly by making the number k of boxes, and their sizes N„, large enough, 
we can approximate any density profile p{x) via p{x) = M p /N p for x = p/k. 
Equation ( [4.9| ) then gives the probability of the profile p(x), and one may 
attempt to obtain an asymptotic form via a saddle point analysis. This 
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approach, which is conceptually straightforward and was followed in the 
symmetric case turns out to be more difficult to implement in the 

asymmetric case, and this is why we here follow the additivity approach 
explained in Section [5|. 



Various more explicit expressions may be extracted from (4.6) and (4.7); 
for example, in the totally asymmetric case (q = 0) the average profile (rj) 
boundaries was computed for all i in | |12[ |. Moreover, for q = and at the 
special point a = (3 = 1 {p a = 1 and pb = 0) in the maximal current phase C, 
finite-size corrections have been computed |2£| for the mean density profile 



(t{) and for the two point function (t^t.,) — (Tj)(r 7 ). In particular it was 



found in [28] that at a point i = Nx (with 1 < i < N) of a system of N sites 



^ = \ + ^70^) + O(N ~ 3l2) (4 ' u) 

which for a system of N = L(b — a) sites becomes for i = L(x — a) 

< T '» - 5 + T7;Tm m "H'X \ + 0{L ' m) (412) 

2 2^ir L l / Z y/(b - a)(b - x)(x - a) 

Also, at this point (a = (3 = 1) of the phase diagram it was shown [pq | 
that for large A r , the variance of the total number M of particles is given by 

(M 2 ) - (M) 2 ~ — (4.13) 
8 

while a Bernoulli distribution at density 1/2 would give twice this variance 
(see the discussion in Section |3.5| ) . 



5 Derivation of the additivity formulae 

In this section we obtain the additivity formulae ([O]) and ( |1.10| ) for the large 
deviation functional in the ASEP. We begin by deriving two formulae valid 
for arbitrary system size N. The first of these, ( 5. ICQ , is obtained directly 
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from the matrix formalism of Section |3| and is valid for a range of parameters 
corresponding to the fan region p a > Pb- By analytic continuation of ( |5. 10Q 
we then obtain a new formula ( 5.23j ), which is valid for all parameter values. 



Finally, we analyze the large N behavior of ( |5.23j ) to obtain (1.6) and Ql.icp . 



For q = 0, the additivity formula ( 5.1Cp takes a much simpler form, 



which is given in [31]. We do not reproduce this formula in the current 
paper because the treatment here of the general case < q < 1 requires 
slightly different notations. 

5.1 Preliminaries 



Let D and E satisfy ( [4.1| )- (|4.3| ) with q < 1. We define the operators d and 
e by 

£> = -J-(l + d), E =-^—(l + e); (5.1) 
1 — q 1 — q 

from ( |4.1D these operators satisfy 

d e — q e d = 1 — q. (5-2) 

We also define eigenvectors \z) and (z\ of d and e, for arbitrary complex z, 
by 

d \z) = z \z) , (5-3) 
1 

(z\e = - (z\ . (5.4) 

If X is a polynomial in the operators D and E (or equivalently d and e), 
then the matrix element (zq\X\z\) / {zq\z\) is a polynomial in 1/zq and z\, 
with positive coefficients whenever the coefficients in X are positive; this is 
easily seen since then X is a polynomial in d and e with positive coefficients, 
and using fl5.2|) one can push all the cf s to the right and all the e's to the 
left, maintaining this positivity. 
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Finally, we define the function (f(z) by 

oo ^ 

<p(z) = £ - z n , (5.5) 

n=0 

where the c n are constructed from the recursion 

c = 1 , c n = (1 - <? n ) c n _i. (5.6) 
One can check easily from fl5.5| ) that 

tp(z) - <p(qz) = zip(z) , (5.7) 

so that 



^ = n i — - • ( 5i 



5.2 Exact additivity formula for \zq\ > \z\\ 

The additivity formula which we prove in this section is that if Xq and X\ 
are arbitrary polynomials in the operators D and E, then for 



N > M (5.9) 



one has 



(zolMiki)^) (51Q) 
cfe (z \X \q n z) ( q n z\ (z\Xi\z{) ( z x 



(z \zi) \z 

1 °° n n 
___ V — 



fil) ~^ °n I 2niz (z \q n z) V z J (z\zi) \z 
where the contour of integration is a circle \z\ = R with 

\zi\ < R < \z \ ■ (5.11) 



Proof: Using ( |5.2| ), any polynomial X in D and E can be written in the 
form 

X = Y,A PtP ,eP' cF. (5.12) 
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As ( 5.1C| ) is linear in Xq and X\, it is sufficient to prove it for X = e p od Po 
and X\ = e Pl d p i, and as (zq\ and \zi) are eigenvectors of e and d, one can 
immediately simplify the problem and limit the discussion to the case 



X = d po 



X x 



(5.13) 



For this choice of Xq and Xi, the right hand side of (5.10) becomes 

oo oo oo n ( 1+mo+po ) mi 



-j^ iA; y^-i 

EE E 



<p(q) 



1 X 

mo u mo+po,mi+pi > 



(5.14) 



once if(q n z/zo) and <p(zx/z) have been replaced by their power series ( |5.5| ) 
and the integration over 2 has been carried out. One can now use (5.5) again 
to evaluate the sum over n, and then simplify the result using the identity 
ip(q 1+m ) = c m (p(q), which follows from ( |5.7| ) by an inductive argument, to 
obtain 



r.h.s. of (5.10) 



00 00 

E E 

mo=0 mi=0 



1 



,mi 
1 A 

m °mo+po,mi+pi) 



00 co mi 
c mo+po z l 



E E 



mo Om +po,"il+Pl 



(5-15) 

m =0ra 1 =0 Lm » Lmi ^0 

So to prove ( 5. ICQ , we just need to show that the left hand side of ( |5,10| ) 
coincides with (|5.15|) . 

We argue by induction on po- If Po = 0, it is easy to see that the left 
hand side of ( 5,10| ) is given by z Pl ip(zi/zo) . Evaluating (5.15) in this case 
leads to the same expression. 

Let us assume, then, that ( 5.1C| ) is valid for all po < P. We want to 
prove that it remains true for po = P + 1. To do so we observe the following 
consequence of (5.2): 



[l-q pi )d p e pi ~ l + q pi d p e pi d. 



(5.16) 
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Using ( EU6[) , the left hand side of ( p0| ) for p = P + 1, i-e., for X = d p+l 
and X\ = e Pl , becomes 



(zol^ew- 1 !^) 
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-,P1 



{z \d p e^d\ Zl ) fzi 
■if — 

\ z o 



(5.17) 



As we have hypothesized that (|5.10 ) is valid for p$ < P, we can replace each 
term of ( 5.17 ) by the corresponding expressions ( |5.15| ), leading to 



l.h.s. of (5.10) = E 



CO oo mi 
Cmo+P z l 



mo =o mi =o Cm ° Cmi z o 



mo 



[(1 - q Pl )5 mo+ p, mi+Pl -i + zi q Pl 5 mo+ p tmi+Pl } 



oo oo 



- E "T~f^" ^m +-P+l,mi+pi-i , (5.18) 

mo =0 mi=0 

which is identical for p = P + 1 to ( |5.15| ) . This completes the derivation of 



5.3 Analytic continuation of ( |5.10|) 

Recall that ( |5loD has b een established for \zq\ > \z\\, with integration con- 



tour a circle \z\ = R with \z±\ < R < \zq\. However, the left hand side of 
this equation is an analytic function defined for all complex zq and z\ except 
at zq = and at the poles of ip{z\/ zq). In this section we make an analytic 
continuation of the right hand side to obtain an integral representation of 
the left side which is valid for any zq, z\ for which the left side is defined. 
In the final representation we will again integrate over a contour \z\ = R, 
but now R will be allowed to take any value for which the contour does not 
pass through singularities of the integrand, that is, for which 



R > 0, 



R ^ q m \zi\, m > 0, R / q~ m \z \, m>0. 



(5.19) 



The expression on the right hand side of ( 5.10j ) must be modified when, 
as one varies zq, z\ and R, the singularities of the integrand — that is, the 
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poles of (p(q n z/zo) and ip{z\/z) — cross the integration contour. When this 
happens, however, the residue theorem tells us how ( |5.10 ) is to be modified: 



one simply includes the residue of the pole on the right hand side of the 
equation, adding or subtracting it according to whether the pole crosses the 
contour from the inside to the outside or vice versa. Using the fact that the 
residue of ip(z) at z = q~ m , m = 0, 1, . . ., is 

f_1 \m+l m(m-l)/2 

<p(q), (5.20) 

Cm 

one then finds the following extension of ( [3.10 ): 



Suppose that R satisfies ( |5.19 ). If R < \z\\ then define k\ by 



q kl+1 \ Zl \ < R < q kl \zi\ , (5.21) 
and if R > \zq\ define fco by 

N < R < W . (5.22) 
qk qk +i y ' 

Then 

(zq\zi) \W 

^ c n \ <p(q) J\z\=R liviz (z \q n z) ^yzo) {z\zi) 

fcl ( i \m.m(m+l)/2 /„n+m„ \ 

{z \X Q \q n+m Zl ) (q m z 1 \X 1 \z 1 ) 



(z \q n + m Zl ) (q m z 1 \z 1 ) 

" l, ~ n (_-|\m^m(m+l)/2 f q n ^~ m Z\ 



W "/ m ra+1/2 /_n+m„ \ 

+ e(l* - g-kol) £ <P r— 3 -) 

(5.23) 



x 

m=0 

(zo|^o|<r"^o} (g- (n+m) z |^iNi> 



(zok" m -2o) (g-("+ m )z kl) 
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5.4 Asymptotics of the additivity formula 

We now ask what happens to the representation ( [5.23D when the system 
size N becomes very large. Throughout this section we take zq and z\ to 
be positive real numbers and Xq and X\ to be polynomials, with positive 
coefficients, in the operators D and E. We will use throughout this section 
three properties of the product 

{zo\X»\ q »z) (z\Xi\») 
nK ' (z \q n z) {z\z x ) y ' 

occurring in the integrands of the representation ( 5.23| ), which follow from 

the discussion of Section 5.1: n n (z) is a polynomial in q n z and 1/z, with 

positive coefficients; for z on the positive real axis, H n (z) is a convex function 

of z; for fixed positive z, TL n (z) decreases as n increases. We will assume 

that, for z on the positive axis, IL n {z) grows exponentially in N. 

Now IIo(z) has a unique minimum at some value z m i n on the positive real 

axis; we will use the representation (|5.23|) with the choice R = z m \ n . Using 

first the fact that the maximum of the magnitude |II n (z)| on the contour 

\z\ = -2min occurs on the real axis, and then the monotonicity of H n (z) in n 

for z positive, we see that each integral occurring in ( |5.23| ) satisfies 

dz (z \X \q n z) (q n z\ (z\X x \ Zl ) ( z x 



\z\=z min 2 ™z (z \q n z) \ z J (z\zi) 



z 



< 



^min) (^min 
(^0 1 Q n Zmm) (^min|^l) 



^ (^ol^ol z m'm) (^miiil^l \ Z X) / c nr\ 

— / i v 1 i \ ) {o.ZOj 
\ZQ\Z m \ n ) \^min|^l/ 

where the first inequality holds up to factors that do not grow exponentially 
with N. On the other hand, when n = the point z m i n will be a saddle 
point for n (z) lying on the contour \z\ = z m - m , and equality (again up to 
factors not growing exponentially with N) will hold in (|5.25 ). 

The bound and an argument similar to that above for the terms 

in ( |5.23 ) arising from the poles, show that the n = terms there always 
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dominate those with n > 0, so that we may neglect the latter. But for 
n = 0, the fact that Ho(z) is convex, with a minimum at z mm , implies that 
the m = terms dominate each sum over m. Thus ( |5,23| ) becomes 



(go^oXi^i) 

(2b | 21 ) 



~ max< 



e(zi 

6 (z m 



(z \X Q \z min ) (Zminl^lkl) 
(^ol^min) (^minl^l) 

(zo\Xo\ Zl ) (zilXijzi) 



2min) 
- 2 ) 



(2 |2l) (2l|Zi) 

(z |Xo|z ) (z |Xi|zi) I 



(5.26) 



(2q|2o) (2o|2l) J 

Now the discussion can be completed by considering successively all the 
possible relative positions of Zq, z% and Zmin- 



Suppose first that Z\ < zq. Then ( 5.26 ) gives 

(z \X \z min ) (^minl^lkl) 



(2d 1 2l) 



(2o|-2min) (2min|^l) 

(z \X \ Zl ) (^JjCiJzi) 
(20 1 21) (21I21) 

(z |X |z ) (2o|^l|2l) 



(2o|2 ) (20^1) 

and one sees that each case in ( 5.27] ) reduces to 



if Z\ < Z min < Zq, 



if z min < zi < zq, (5.27) 



if z\ < zq < z n 



(zolXoXxlzx) 



~ mm 



(zq\X \z) (z\Xi\ Zl ) 



(zo\z\) Z!<z<z (zq\z) 

On the other hand, if zq < z\, then ( 5.26Q gives 

(zq\X \z) (z|Xi|*i) 



{zolXpXM 
(20I21) 



max 

2=20, 21 



(zo\z} (z\zi) 

(£q|Xq|£i) (z^X^Zj) 
{z \zi} (Zl\z!) 

(z \X \zo) {z \Xi\zi) 



(2o|2 ) {z \zi) 
and each case in ( fgp reduces to 



(2I21) 

if ^0 < 2 < Zi, 
if 2 min < 20 < 21, 

if Z < 2i < Z min , 



(5.28) 



(5.29) 



(zo|jr *l|*i) (2 |X |z) <*|Xi|si} 

~ max 



(zo 1 21) 



2 = 20,2! (Zo|2) (2|^l) 



(5.30) 
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5.5 Derivation of ( |1.6| ) and ( |1.10| ) 



We finally want to obtain the fundamental additivity relations for 7i, ( |1.6| ) 
and ( |1.10| ), from (|5.28|) and ( |5.30 ). The first step is to relate the densities 
Pa, Pb to the parameters zq,z\. Using ( |5Ji| ) and (^4|), one can easily 

establish that if 

z z 1 



then 



Pa : 

(zq|XqXi|2:i) 

(zq\zi) 



Pb : . 

1 + Z\ 

(WIXqX^V) 
(W\V) 



(5.31) 
(5.32) 



where (W\ and \ V) are defined by (|4.2|) and (|4.3| ). 

Now let us consider a given profile /o(x) defined for a < x < b, and for 
fixed c with a < c < b denote by Xq the sum over all the products of -D's 
and E's consistent with the left part of this profile over the first L(c — a) 
sites, and by X\ the same quantity for the right part of the profile over the 
last L(b — c) sites. We define TL by 

\L(b-a) { z o\ X Q X l\z\) 



(i - qY 



(zq\zi) 



exp [-LH[ ajb ]({p(x)};p a ,p b )] . (5.33) 



Then we obtain immediately fll.6j ) and Ql.lCD from ( |5.28| ) and ( |5.30| ). More- 
over from (1.2), (4.4), Q5.32 ), and ( 5.33] ) we see that the constant K(p a ,Pb) 
which appears in ( |1.4j ) is given by 

(5.34) 



(1 - q) L {b - a)(W\{D + E) L {b - a)\V) ~ e - L (b-a)K{p a , Pb ) 



Writing (ET28D and (|53o|) for X = (£> + £) L(c ~ a) and X x = (£> + E) L ( b ~ c ), 
we see that K{p a ,pb) should satisfy 

(b - a)K(p a , p b ) = sup [(c- a)K(p a ,p c ) + (b- c)K(p c ,p b )} , (5.35) 

Pi)<Pc<Pa 



if p a > Pb, and 
(6- a)K(pa,Pb) 



^ in ™- „ [(c-a)K(p a ,p c ) + (fo-c)^^,^)] » 

" Pa Cil pfe 

(5.36) 
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if p a < pb- Now when p a = p the matrices D and E commute and may be 
realized as the scalars D = (1 — + zq), E = (1 — + z^ 1 ), so 

that from ( 5.34 ) and ( 5.31| ) we have 



K(p a ,Pa) 



log ( ZQ + 2 + — ) = l0g[p o (l - p a )} . 



(5.37) 



From ( |5,35| ) and ( |5.37| ) one finds, by repeated subdivision of the interval 
[a, b], that if p a > Pb then 

(b - a)K(p a , Pb ) = sup f b dx log[/(x)(l - f(x))} , (5.38) 
/ Ja 

where the supremum is over nonincreasing functions f(x) with f{a) = p a 
and f(b) = pb, and from ( |5.38; ) one obtains (|1.5|). Similarly, if p a < Pb one 
obtains from ( |5.35| ) and ( |5.37| ) that 

r rv rb 
(b - a)K(p a , p b ) = inf <^ / dx log[p a (l - p a )} + / dx \og[p b {l - p h ) 

a<y<b {Ja Jy 

(5.39) 

and (lL9|) follows. 



6 Large deviations versus typical fluctuations 

Expressions (1/7) and ( l.llj ) enable us to calculate the large deviation func- 
tion for an arbitrary density profile p(x). As we have already noted and 
will show in this section, the large deviation functional, which describes 
macroscopic (order N) deviations from the typical profile p, is in region C, 
where 

Pa > ^ > Pb , (6-1) 



not simply related to the fluctuations (order vN) around p. We demonstrate 
this by computing explicitly, for q = 0, the probability of seeing a given 
global density r in a window c < x < d of our system, with no other 
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constraint in the system, both for fixed r with r ^ p and for r — p of order 

i/Vn. 

We divide our system of N = L(b — a) sites into three boxes, of N± = 
L(c — a) sites, N2 = L(d — c) sites and N3 = N(b — d) sites, corresponding 
to macroscopic intervals [a, c], [c,d], and [d,b], and compute the probability 
that the density is r in the middle box, i.e., the probability P^My) that the 
total number M2 of particles in the middle box is 

M 2 = L(d-c)r, (6.2) 

with no constraint imposed in the two other boxes. 

6.1 Large deviation 

Corresponding to the above constraint there will be an optimal profile p(x) 
for which, with T(r) = F[ a ,b](p( x )] Pa, Pb), one has 

P(M 2 )~exp[-L^]. (6.3) 

Since p a > pb, one can use the additivity formula (|1.6|) 

W[a,6]OKaO;Pa,P&) = SU P {^[a,c](p(x);p a ,Pc) 

Pb<Pd<Pc<pa 

+ M[c,d](p(x);Pc,Pd) + {H[ d)b ](p(x); p d , p b )} • (6.4) 

As there is no constraint on the profile in the two side boxes, one has 
F[a,c](p&); Pa, Pc) = F[d,b](p{x); Pd, Pb) = 0, so that from (|l~4|) , 

Ti[a,c\{p{x); Pa, Pc) = (c- a)K(p a ,p c ), (6.5) 
7-L[b,d](p(x); Pd, Pb) = (b- d)K(p d ,p b ). (6.6) 

Moreover, since p c > pd we know from Sections |3.5| and |3.6| that within 
the central box [c,d] the optimal profile p(x), which corresponds to a fixed 
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number of particles there, is flat. From (1.4) and (3.12) we then have 

(c - d)[rlog(r(l - p c ) + (1 - rlog((l - r)p c )], 
(c-d)[2rlogr + 2(l - rlog(l - r)}, (6.7) 
k (c- d)[rlog(r(l - + (1 - rlog((l - r)p d )], 



when r < 1 — p c , 1 — p c < r < 1 — p d , and 1 — < r, respectively. Now 
from ( |6.4| ) one must choose p c and to maximize the sum of (|6.5|), (|6.( 
and ( |6.7| ). The result depends on the sign of r — 1/2. 

Case 1: r < 1/2. In this case the optimizing values of p c and pd are 



id — a) — rid — c) 
Pc = mm< < d+c-2a > Pa 



Pd 



(6.8) 



The corresponding LDF is 



F(r) = (d-c)[rlog(4r(l-p c )) + (l-r)log(4(l-r)p c )] 

+ (c-a)log(4 Pc (l- /0c )), (6.9) 



which for r close to 1/2 becomes 



4 {i -c)( d -a)f _iy / 
w d + c-2a\ 2 y V 2 



(6.10) 



We also find from Section 3.2 that the optimal profile p(x) satisfies p(x) = p c 
for a < x < c and p(x) = 1/2 for d < x < b. 

Case 2: r > 1/2. In this case, a similar calculation leads to 



1 



Pc 



p d = max 



- c) — r(d — c) 
2b-d-c 



<Pb 



(6.11) 



and 



F{r) = (d-c)[rlog(4r(l-p d )) + (l-r)log(4(l-r)p d )] 

+ {b - d)\og(Ap d {l - p d )); (6.12) 
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and for r close to 1/2, ( |6.12j ) gives 

We see from ( |Sl0D and flOfl that in general the limiting value of the 
second derivative of T{r) at r = 1/2 depends on the sign of r — 1/2, so that 
J-{r) is in general nonanalytic at r = 1/2. Note, however, that when the 
overall density for the system is specified, i.e., when c = a and d = b, J-(r) 
is analytic. 



6.2 Fluctuations 



Formulae ( |6.9[ ) and ( 6.12| ), with ( |6.8[ ) and ( 6.11| ), give us the leading behavior 
of P{M-z) for large deviations, i.e. for r — 1/2 of order 1. Let us now discuss 
the small fluctuations, i.e., the regime in which r - 1/2 = 0(L~ 1/2 ), so that 
LT is of order one. As ( 6.10| ) and (6.13) do not coincide, one expects that 
typical fluctuations around the optimal profile p(x) = 1/2 will be anomalous, 
i.e, non-Gaussian. Let us define the fluctuation fi in the number M2 of 
particles in the central box by 

(d - c)L 



Mo 



fiV L. 



(6.14) 



In the next subsection we will derive, for q = 0, the following probability 
density p(n) = P(M<2)\fL for the random variable \x: 

8 (b - af' 2 



p(m) 



[ir(c-a)(b-d)] 3 / 2 (d 

00 />oo 

dx I dy xy exp 
(I Jo 

2 fJ 2 + {fi + x - y) 



x 



(6.15) 



xr 



2 1 



exp 



c — a 

2' 



6-d 
exp 



(/i + x) 2 + (ji - yf 
d — c 



We see that indeed \i has a non-Gaussian distribution. Note, however, that 



the V L scaling in ( 6.14 ) is that appropriate for "normal" fluctuations. 
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For large positive //, ( |6.15| ) becomes 

4 {b-af/ 2 (b-c){b-d) 



1-iTT (d - c) 1 /2(26 - c - dfl 2 (c - afl 2 

4(6 - c) 



exp 



and for large negative fj,, 

p{p) : 



(d - c)(2b - d - c) 
{b-afl 2 {d-a){c-a) 



(6.16) 



(6.17) 



(d - C y/ 2 (d + c - 2a) 3 / 2 (6 - d) 3 / 2 
T 4(d - a) 

6XP [~(d-c)(d + c-2a)' 

We see that the large /U asymptotics ( 6.16D and ( 6.17| ) match with ( |6.10| ) 

and ( |6.13| ) when \n\ S> 1 and \r — 1/2 1 <C 1, with fi = (r — l/2)y/L, i.e., 

that ( |6.15| ) interpolates between the large deviation regions r > 1/2 and 

r < 1/2. A similar relationship has been found between the distributions of 

large deviations and of typical fluctuations of the current in the ASEP on a 



ring [32]. 



Also, from (6.15) one can compute the average (/z) of fi: 

y/(d - a){b -d)- y/(c - a){b - c) 



For a system of N = L(b — a) sites, one sees from (4.12) that 



L{d~a) 
i=L(c—a) 



1\ _ y/L f d b + a-2x 

{Ti} ~2)~2^J C y/(b-a)(b-x)(x-a) 



dx, 



(6.18) 



(6.19) 



in agreement with ( |6.18| ), since (M2) — L(d — c)/2 = (/^)VX. 



6.3 Derivation of (|6.15|) 

One way to derive fl6.15| ) is to use the following explicit representation 



of flOD, valid for 9 = 0: 

00 

D = + |ra)(n+ 1|, 



(6.20) 



n=l 
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E = |ra)(n| + |n + 



(6.21) 



n=l 



where the vectors |1), |2), ...|n)... form an orthonormal basis of an infinite 
dimensional space (with (n|m) = 5 n ^ m ). Within this basis, the vectors \V) 
and (W\ are given (see Q, Q, and (|Ll|)) by 



1*0 = £ 

n=l 

oo 

w = E 



1-/3 

/? 

1 - a 

a 



|-> = E i 

n=l 

oo 

H = E 



Pb 



Pb 

1 - Pa 
Pa 



n=l v x n=l 

and one can show (for example by recursion) that 

(2AQ! 



n) 



(6.22) 
(6.23) 



and that [j33| l 
(p|Xjv,Mb') 



(2iV)! 



(iV+p-p')! (iV+p'-p)! (iV+p + p')! {N -pf -p)\ 

(6.24) 



(AT! 



n 2 



(M)! (JV — M)! (M + p- p')\ (N -M-p + pty. 
{N\) 2 



(M + p)\ (N - M - p)\ (M - p?)\ (N -M + pf)\ ^ 6 ' 25 ^ 
where Xj^,m 1S the sum over all the configurations of ./V sites with M occu- 
pied particles. The probability P(M 2 ) that the number of particles is M 2 in 
the central box is given by 

(W\(D + E) N i X N2 , M2 (D + E) N *\V) 



(W\(D + E) N i+ N 2+ N a\V) 



P(M 2 

Let us first analyze the denominator of fl6.26j ): 

1 - Pa\^ ( Pb 



(6.26) 



{W\{D + E) N \V)=Y J E 

Pl=lp2 = l 



Pa 



I- Pb 



( Pl \(D + E) N \p 2 ). 



(6.27) 

For large N, this sum is dominated by p\ and p 2 of order 1, so one can use 
an approximation of ( 6.24 ) valid for p and p' of order or less, 



(p\(D+E) N \p" 



cxp 



'\2 



N 



exp 



■jp + P 
N 



l\2 



, (6.28) 
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4^ (l- Pa )Pa (l-p b )Pb 



and one gets for large N, 

One can write the numerator of fl6.26j ) as 

(W\ (D + E) Nl X N2 , M2 (D + E) N *\V) =EEEEW« 



(6.29) 



PI P2 P3 P4 
A^ 3 I 



( P1 \(D + E) 1 ^ \p 2 ) ( P2 \X N2 , M2 \p 3 ) (p 3 \(D + E) N * \ P4 ) (p 4 \V). (6.30) 

When N\, N 2 ,N 3 are large and of order L and when the difference M2—N2/2 
is of order y/L as in (1P3) , these sums are dominated by pi and ^4 of order 
1 and P2 and j»3 of order y/L. If one writes 



one gets that 

{P2\X N2 ,m 2 \pz) 



p 2 = xyfZ, 



P3 = y 



Vl, 



(6.31) 



vriV 2 
— exp 



cxp 



2L 

- — {p 2 + {p + x - yf) 

2 ^-y) 2 + ^ + x?) 



,(6.32) 



and the numerator becomes, after summing over p\ and p 4 and replacing 
the sums over P2 and p 3 by integrals 



(W\(D + E) N > X N2 ,m 2 (D + E) Ns \V) 

^Ni+N 2 +N 3 +2 x 2 



(6.33) 



TT 2 L 2 [{c-a){b-d)f/ 2 {d 
(l-p a )p a (l-p b )p b 



(2p a -iy (2^-1) 



dx 



dy x y exp 



y 



cxp 



p 2 + (p + x- y) 
d — c 



21 



cxp 



c — a b — d 
{p + x) 2 + {p - y) 2 
d — c 



which reduces to ( 6.15 ) after dividing by Q6.29| ). 

Remark: One can also recover from ( |6.26D the expressions ( |6.9D and 
( |6.12| ) by allowing deviations in M2 of order L. 
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7 Conclusion 



The main results of the present work are the exact expressions ( |1.7| ) and 
(PH) for the LDF F({p}) for the SNS of the open ASEP in one dimension 
and the simple additivity formulae ( |l.6[) and ( 1.10| ) that they satisfy. 



This J-, like the one we found for the symmetric case in [16, |lfif[ , is a 
non-local functional of the density profile {p(x)}. We expect non-locality to 
be a general feature of such functionals for non-equilibrium systems. 

Our expressions of the LDF, which take different forms in the fan region 
Pa > Pb (where the reservoirs and the bulk asymmetry cooperate) and in the 
shock region p a < pi, (where they act in opposite directions), reflect several 
qualitative differences: 

In the fan region, p a > pi,, the probability of macroscopic deviations from 
the typical density profile is reduced compared with that in an equilibrium 
system with the same typical profile (see (|3.9| )); this was also true in the 



symmetric case [16]. Another surprising feature of the fan region, at least 
in the maximal current phase C, is that the fluctuations of the density 
profile cannot be calculated from the LDF, and that these fluctuations are 
in general not Gaussian (see section 6). We have no heuristic explanation 
of this behavior. 

In the shock region of the phase diagram, p a < p^, T is not convex in p 
(see section 3.3). Moreover, in this region the probability of macroscopic de- 
viations from typical behavior is increased rather than reduced; see ( |3.11 ). 
This enhancement of deviations appears similar to known behavior Q of 
fluctuations (the behavior of macroscopic deviations is not known) in a slab 
of fluid in contact at the top with a heat reservoir at temperature T a , and at 
the bottom with another reservoir at temperature T&: the Rayleigh-Benard 
system Q . In this system the force of gravity causes the SNS to undergo dy- 
namic phase transitions when (2], — T a ) is "sufficiently" large, corresponding 
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to different spatial patterns of heat and mass flow. This transition is pre- 
ceded, as Tf, — T a is increased, by enhanced fluctuations even for very small 
differences between the two temperatures (for which the system is stable) 
as long as T& > T a . 

It is natural to expect that the probabilities of typical fluctuations and 
of large deviations are either both enhanced or both reduced, in comparison 
with the equilibrium system having the same p, but this is not known to be 
true in general. In fact, our work here shows that small (typical) fluctuations 
cannot in general be computed from the LDF. 



One can note that our expressions for the LDF (1.7) and (1.11) do not 
depend on the asymmetry parameter < q < 1. These expressions, how- 
ever, are not valid at q = 1, and they do not reproduce the symmetric 
exclusion result: the limits q — ► 1 and N — ► oo do not commute. It would 
be interesting to analyze the case of large N with 1 — q = 0(A^ _1 ), in order 
to interpolate between the symmetric case and the asymmetric case. 

It would also be desirable to have a physical understanding of our additiv- 



ity formulae ( p..6| ) and ( 1.10 ), and to see how our results could be generalized 
to more complicated non-equilibrium steady states. 

Lastly, from our knowledge of the LDF in the SNS steady state, one 
could try to determine how a given (unlikely) profile was produced dynam- 
ically out of the nonequilibrium steady state. This has been done for the 
symmetric exclusion process on a circle in [[34|1, and recently for the open 
system in p5| , |36|| , where the LDF was given in terms of a time integral over 
a trajectory taking the system, via a "reversed dynamics," from a typical 
SNS configuration to the profile p(x). A dynamic LDF for the ASEP on the 
circle was studied in [37]. 
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A The concave envelope construction 



In this appendix we justify the construction (|3.5 ,3.6) of the optimizing func- 
tion F p for the supremum in ( |1.7| ) . Recall that we are given a density profile 
p{x) which is defined for a < x < b and satisfies < p{x) < 1 for all 
x, and reservoir densities p a and pb which satisfy 1 > p a > pb > 0. Let 
H p {x) = (I — p{y)) dy, so that G p is the concave envelope of H p , and re- 
call that F p is obtained by cutting off G' p (x) at p a and pb (see ( |3.6|) ). Then 
(1.7) may be written as 

•^[0,6] ({p(^)>; Pa, Pb) = ~(b- a)K(p a ,p b ) 
+ f dx [p(x) log p(x) + (1 - p(x)) log(l - p(x))} + sup B p (F) , (A.l) 

Ja F(x) 

where 

F(x) 



B P {F) 



dx 



1 - p{x))\0£ 



+ log(l - F(x)) 



(A.2) 



1 - F(x) 

and the supremum in (|A.1[ ) is over monotone nonincreasing functions F(x) 
satisfying F(a) = p a , F(b) = pb- Now since H p {x) < G p (x) for all x and 
log[-F(x)/l — F{x)} is decreasing, we obtain by integration by parts, noting 
that H{a) = G(a) and H(b) = G(b), that for any F(x) in this class, 



dx[(l-p(x))-G'(x))] log 



F(x) 
1 - F(x) 
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[H p (x))-G p (x)]d 



log 



F(x) 
1 - F(x) 



< 0, (A.3) 



so that 



B P {F) < f 

J a 



dx 



G' p {x)\0g - F{ l)^ -Mooil - FUr)) 



(A.4) 



1 - F{x) 

Since the integrand in ( A.4 ) is pointwise concave in F(x), with a maximum 
at F(x) = G' p (x), the integral is, for the functions F(x) satisfying p a > 
F(x) > p)b, bounded above by its value at F{x) = F„(x). Thus 



B P (F) < J" dx G' p {x) log - i + log (i _ j? 



c/.T 



d -/>'■'•)) Fp i x L +iog(i-w) 



^p(^) ; 



(A.5) 



here the first equality is obtained by noting that if (c, d) is a maximal interval 
on which G p (x) ^ H p {x) then (i) G' p (x) and hence F p (x) are constant on 
this interval and (ii) f^G' p (x)dx = J^(l — p{x))dx. Equation ( |A.5| ) shows 
that the supremum in ( |A.1| ) is achieved by F(x) = F p {x). 



B Optimal profile under a constraint 

Let p(x) be an optimal system profile for a fixed mean density r, that is, a 
profile which minimizes F{ a fi} ({p( x )}l Pa, Pb) under the constraint 



p(x) dx = r(b — a). 



(B.l) 



In this section we show that p(x) is uniquely determined, and derive its 
form. 

Case 1. p a > Pb- Let G p and F p be defined by (|3.5[) and (|3.6|) , so that 
F[ a ,b]({p(x)};Pa,Pb) = h(p(x),F p (x);p)dx, (B.2) 

J a 
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where h(r,f;p) is defined in ( |3.2| ) and p is determined from p a ,Pb as in 
Figure 1. 

We first show that p(x) must be monotone nondecreasing, that is, that 
p(x) = 1 — G' p (x) (almost everywhere). For otherwise there will exist some 
interval [c, d] C [a, b] satisfying 

G p (c) = f (1 - p(y)) dy, G p (d) = f\l - p{y)) dy, (B.3) 

J a J a 

and 

rx 

G p (x) > / (1 - p(y)) dy, for c < x < d. (B.4) 

J a 

Then for c < x < d, G' p {x) = 1 — A, where 



1 



A = — ^ / p(y)dy, (B.5) 



and if p* is defined by 

p*(x) 



A, if c < x < d, 
p(x), otherwise, 



(B.6) 



then p* satisfies (BA) and from the strict convexity of h(p,F;p) in p and 
the fact that F p is constant on [c, d] (with value A, p a , or pt,) it follows that 

(c-d)h(X,F p ;p) < J h(p(x),F p (x);p)dx (B.7) 

and hence that T[ a ^{{p* (x)}\ p a , p b ) < T[ a ^{{p{x)}\ p a , pb), contradicting 
the definition of p(x). 

Since p(x) = 1 — G' p (x) we have from (|3.6|) that F p (x) is a local function 
of p(x) taking value p a if 1 — p{x) > p a , p^ if 1 — /)(x) < and 1 — p{x) 
otherwise. Then one can check that the integrand h(p(x), F p (x); p) in (|B.2| ) is 
pointwise convex in p(x). Thus since p{x) satisfies (|B.1| ), (b—a)h(r, 1— r, p) < 
^\a,b}{{p{ x )}'^ Pa-, Pb), with equality if and only if p(x) is the constant function 
p(x) = r. 
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Case 2. p a < p b . Let y r be the minimizing value in ( |1.11 ) corresponding 



to the optimal profile p(x); it is clear that p(x) must be constant on the 
intervals [a,y r ] and [y r ,b], so that 

F[a,b]({p( X )}'' P°» Pb) ( B - 8 ) 
= (y r - a)h(r a ,p a ; p) + (b - y r )h(r b ,p b ; p) 

for some values r a ,r b satisfying 

(y r - a)r a + (b - y r )r b = (b - a)r. (B.9) 



Minimizing ( |B.8[) over a < y r < b and ( [B.9D leads to y r = b, r a = r if 
r < 1 — p b and y T = a, r b = r if r > 1 — p a , while if 1 — p b < r < 1 — p a , J/r 
is given by ( ^.20[) and r a = 1 - p b , r b = 1 - p a . 
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